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ABSTRACT 

During two weeks of nearly continuous optical photometry of the young open cluster 
^ . NGC 2244 obtained by the MOST satelHte, we discovered two new SPB stars, GSC 00154- 

; 00785 and GSC 00154-01871. We present frequency analyses of the MOST light curves of 

, these stars, which reveal two oscillation frequencies (0.61 and 0.71 c/d) in GSC 00154-00785 

and two (0.40 and 0.51 c/d) in GSC 00154-01871. These frequency ranges are consistent with 
'nI" ' g-modes of ^ < 2 excited in models of main-sequence or pre-main-sequence (PMS) stars 

. of masses 4.5 - 5 Mq and solar composition {X,Z) - (0.7,0.02). Published proper motion 

measurements and radial velocities are insufficient to establish unambiguously cluster mem- 
bership for these two stars. However, the PMS models which fit best their eigenspectra have 
ages consistent with NGC 2244. If cluster membership can be confirmed, these would be the 
first known PMS SPB stars, and would open a new window on testing asteroseismically the 
interior structures of PMS stars. 

\ Key words: stars: early-type - stars: pre-main sequence - stars: individuakGSC 00154- 

■ 00785 and GSC 00154-01871-stars: oscillations - open clusters and associations: individual: 

NGC 2244 



1 INTRODUCTION 

Slowly pulsating B (SPB) stars pulsate in high-order, low-degree 
gravity modes with periods of the order of a few days, show- 
ing dense frequency spectra, low amplitudes and multi-periodicity 
with beat period s of months and even years ( Wael kens et al.|[l"987l : 
IWaelkensllT99ll) . The effective temperature range of known SPB 
stars is 10,000 - 20,000 K; the mass range is 3 - 8 solar masses. SPB 
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pulsations are driven by t he <c-mechanism of t he Fe opacity bump 
at T ^ 1.5 X lO^K (Gau tschv & Saidll993l : iDziembowski et ai.l 
1 19931) . Being trapped deep inside the star, SPB g- modes are ver y 
promising for asteroseismology of massive stars toe Catll2007h . 
Currently, there are at least 50 confirmed SPB stars and 65 can- 
didate SPB stars in o ur Galaxy, 60 in the LMC and 1 1 in the SMC 
l lKolaczkowski et al.ll2 006). 

Since SPB stars have periods of a few days, it can be diffi- 
cult to characterise their eigenspectra accurately with data from 
ground-based observatories due to aliasing caused by daily gaps 
in the time series. Space-based photometry, with long, nearly unin- 
terrupted time coverage, is essential to perform the best frequency 
analyses for asteroseismology. 

The MOST (Microvariability & Oscillations of STars) mis- 
sion, Canada's first space telescope, obtains highly precise opti- 
cal photometry of bright stars (V < 10") nearly continuously over 
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time spans up to 60 days jWalker et all2003h . Within two years of 
MOST's launch in 2003, the satellite's capabilities were expanded 
to include photometry of the Guide Stars for each Primary Science 
Target field. Many variable stars have been discovered in this large, 
nearly random. Guide Star sample (9 < V < 6), spread across the 
HR Diagram. 

In January 2008, MOST observed the very young open cluster 
NGC 2244. Embedded in the Rosette Nebula, NGC 2244 is located 
in the vicinity of the Galactic anti-center (/ = 206.306 b = -2.072) 
at a dis tance estimated to be in the range of 1.4 < D[kpc] < 
1.7 ( see 'Ogura&Ishida"l98l'; 'Perez et al."l987|; iHensberge etal] 
12000 : Park & Sung 2002; Bonatto & Bica 200%. 

The a ge of the cluster was first estimated to be 2 - 4 Myr (see 
IChen et aU2007. and references ther ein); a more recent esti mate al- 
lows a broader range of 0.2 - 6 Myr teonatto & Bicall2009l) . While 
many of the cluster stars appear to be younger than 3 Myr, about 
45% of the cluster members may be older than 3 Myr, indica ting 
an earlier epoch of star formation teergofer & Christianll2002l) . Of 
the younger portion of the sample, 5% are estimated to be signifi- 
cantly younger than 3 Myr, suggesting star formation is ongoing in 
NGC 2244. 

GSC 00154-00785 and GSC 00154-01871 were used as Guide 
Stars for the NGC 2244 field and were found to be variable. The 
data analyses presented in this paper show the stars to be SPB stars 
similar to HP 16 3830, the first SPB star discovered by MOST 
jAerts etal.l l200^ In SectH we give a short description of the 
photometry, and in Sect[3l our analysis techniques and results. 
Sect|4]deals with evidence for the cluster membership of both stars. 
Section [5] describes the theoretical models which best fit the ob- 
served pulsations, and we present discussion and our conclusions 
in Sect|6] 



2 OBSERVATION AND DATA REDUCTION 

The MOST satellite monitored GSC 00154-00785 (B9V, V= 10.93, 
a = 6''33"'35.7', S = 5°16'41") and GSC 00154-01871 (B5V, 
V=11.30, a = 6''33"'43.3', S = 4°55'54") as guide stars for sci- 
ence observations of the NGC 2244 cluster field. Fig.[T]and[2]show 
the light curves of GSC 00154-00785 and GSC 00154-01871, re- 
spectively. Both light curves clearly show several simultaneous os- 
cillations with periods near 2 days. 

For the NGC 2244 observations 30 subsequent exposures were 
co-added ('stacked') onboard the satellite with a respective expo- 
sure time of 2.03 s. This results in a total sampling period, i.e. 
the time between consecutive photometric measurements, of 60.9 s. 
GSC 00154-00785 was guide star #35 and GSC 00154-01871 was 
guide star # 36 during this observation run. 

Stray light due to scattered Earthshine reaches the focal plane 
of the MOST Science CCD, and this background is modulated with 
the orbital fr equency of the satel lite (14.2 c/d) and some of its har- 
monics. (See Reegen et al. '2006' for a discussion of this effect.) The 
Guide Star photometric data were folded in phase at the period of 
the MOST orbit (about 101 min) to treat the orbital modulation of 
the stray light background. Data points most strongly affected by 
stray light and other outliers were rejected from the light curve. 
To account for the changing nature of the Earth's albedo beneath 
the moving satellite, this process was performed with consecutive 
segments of the light curve, each four MOST orbits long. 

We draw attention to a decrease in relative flux of about 0.02 
mag apparent in all the stars in the NGC 2244 field during JD- 
JD2000 = 2936 - 2940 (and noticeable in Figures □ and O. Al- 
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Figure 1. Light curve of GSC 00154-00785 obtained by MOST photometry 
for 21 days after subtracting the mean magnitude. The red solid line shows 
a boxcar running average with a smoothing length of 50 s. 
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Figure 2. Same as Fig.[T]but for GSC 00154-01871. The red solid line 
shows a boxcar running average with a smoothing length of 200 s. 

though this artefact does not affect our pulsational frequency anal- 
yses below, we do conservatively neglect any frequencies in the 
time series with / < 0.25 c/d, since we cannot be sure they are 
intrinsic to a given star. 



3 DATA ANALYSIS 
3.1 Time series analysis 

For o ur frequency search, we employed PeriodO'^ jLenz et akl 
I2OO5I1 . Period04 is a code specifically designed to extract multiple 
periodic signals in astronomical time series through simultaneous 
least squares fitting. 

The assessment of the signal-to-noise ratio (S/N) for a given 
amplitude in general is a complex task. With Period04, a mean 
amplitude is determined in an interval centered on the frequency 
of interest which is then compared to the amplitude of the peak 
in question. Obviously, the result depends on the chosen interval 

' http://www.univie.ac.at/tops/Period04/ 
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GSC 00154-00785 






/ [c/d] 


<T amp [mmag] 


h 


0.608 


± 0.048 


38.2 


fi 


0.705 


± 0.048 


32.0 


h 


(/i + fi) 


± 0.048 


6.9 




GSC 00154-01871 




fi 


0.407 


± 0.048 


17.6 


fi 


0.514 


± 0.048 


8.20 



Table 1. Frequencies and combination frequencies of GSC 00154-00785 
and GSC 00154-01871 with errors estimated according to o" = T"', and 
amplitudes (amp). 



which may co ntain non-noise signal due to intrinsic frequencies, 
side-lobes, etc. Vau ghanI introduced a method (used by e.g. 

[Oruberetal. 2011) with which one can estimate whether quasi- 
periodic-pulsations (QPPs) are indeed significant. We applied this 
method to check which amplitudes exceed the 3.8 tr limit com- 
monly used in asteroseismology and find 3 and 2 frequencies pass- 
ing this requirement for GSC 00154-00785 and GSC 00154-01871, 
respectively. 

To determine the uncertainties, we use the well known 
Rayleigh criterion {cTf = T"', with T being the data set length), 
which is a cons ervative estima t e of t he frequency error as was 
demonstrated bv lKallinger et al.l j2008h . 

Table [T] lists the oscillation frequencies and uncertainties, 
as well as the associated amplitudes that we obtained using the 
methodology outlined above. Figures[3]and|4]show the correspond- 
ing amplitude and power density spectra. The number of identified 
frequencies i s small compared to other know n SPB stars, as e.g . 
HD 163830 jAerts et al.ll2006l). H P 163868 ( iWalker et al.ll200g) . 
and HD 163899 jSaio et al.ll2006l) . This deficit may be due to a 
larger noise level of our observations in the frequency domain due 
to a data set, which is approximately half as long as that for the 3 
mentioned SPB stars. Additionally, the resulting lower frequency 
resolution reduces the capability of identifying closely spaced fre- 
quencies. 

iDe Cat! ( |2007|. and references therein) argues that amplitude 
spectra of SPB stars with only a few dominant modes can be 
explained partially or even entirely by stel lar rotation, a con- 
cept which was elaborated more recently by i Balona et al.l bOllh 
for yDor stars. If future observations of GSC 00154-00785 and 
GSC 00154-01871 do not unravel additional low-amplitude fre- 
quencies, they would indeed be candidates for such an explanation. 
However, if one considers the observed v sin ( of about 8 km/s (see 
Section [321l, one finds that GSC 00154-00785 and GSC 00154- 
01871 would need to have an inclination angle, of about 4 to 5 
degree in order to explain the observed frequencies due to rotation. 
Assuming a random distribution of ; the probability of finding a 
star with ; < 5 is P = 1 - cos(5) which translates to P » 0.4%. 
At this point we cannot rule out this explanation for our observed 
frequency spectra. But we stress that it is rather unlikely that rota- 
tional modulation is the cause for the observed periodicities. In any 
case, SPB stars with few pulsation frequencies, or few frequencies 
dominating significantly the frequency spectrum pose a challenge 
to theory and provide an interesting observational fact in itself. 
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Figure 3. Upper panel: Region of interest of the amplitude spectrum of 
GSC 00154-00785 (blue) with the identified frequencies (black). Lower 
panel: Amplitude spectrum in log - log space. Shown are the broken power- 
law fit to the data (green solid line) and the 3.8 cr limit (dashed line) 



3.2 Atmospheric parameters. 

For b oth stars Stromgren uvbyf} ph otometry is available i Handled 
I2OIII) . Applying the calibrations bv lNapiwotzki et al] i ll 9931) to the 
data provide the atmospheric parameters which are listed in Table 
2. GSC 00154-01871 is a spectroscopic binary (see lHuang & Giei 
|2006.) , therefore the Hj8 measurement is highly affected by its com- 
panion and thus, not useful to determine log g. 

We also analysed high-resolution (R = 48 000) spectra of the 
two stars obtained during 6-8 December 2008 with the FEROS 
spectrogra ph, installed on the 2.2 m ESO telescope at La Silla, 
Chile (Kaufer et alj|2000l) . FEROS is a cross-dispersed, fibre-fed 
echelle spectrograph. Typical signal-to-noise ratios of the spectra 
(covering wavelengths from 3800 to 9100 A) range from 100 to 
150. From t he FEROS spectra of these stars, provided to us by 
iNiemczural Jpriv. comm.h . we obtain similar parameters to those 
Usted in Table 2, and we confirm that GSC 00154-01871 is indeed 
a spectroscopic binary. From the spectral analysis we also report a 
V sin ( of about 8 km/s for both stars. The corresponding results of 
the spectroscopic observations will be the subject of a forthcoming 
paper (Niemczura et al., in preparation). 

The average atmospheric parameters (from Stromgren pho- 
tometry and from spectroscopy), presented in boldface in Table 2, 
will be used in the remainder of this paper. 
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Table 2. Atmospheric parameters for both GSC 00154-00785 and GSC 00154- 
01871 determined from Stromgren photometry and spectroscopy. 



Star Tetf [K] 


logg[cm/s-] 


ref. 




4.0 ±0.2 
4.3 ± 0.2 


Niemczura (priv. coram.) 
Handler (IQW) 


14600 + 330 


4.15 + 0.15 


average 



GSC 00154-01871 



13800 : 
13605 : 
14300 : 
13900 : 



300 3.8± t.b.d 

59 3.80 ±0.01 

600 4.4 ± 0.2 

360 3.8 + 0.15 



Niemczura (priv. commj 
Huang & Gies (2006) 
Handler ( 2011) 

average" 



for the average of log g the value obtained from photometry was excluded 
because H/3 is strongly affected by the companion. 



Table 3. Proper motion and radial velocity of t he NGC 2244 open clust er an d the two target stars: G SC 00154-00785 and GSC 00154-01871. All values, 
besides the stellar radial velocity, are taken from lKharchenko et alj j2004l) and lKharchenko et alj i2003) . The radial velocity of GSC 00154-01871 is missing 
because the star is a spectroscopic binary and we do not have enough spectra to measure a reliable mean radial velocity. 
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NGC 2244 


06 31 55 


+04 56 30 


-1.42 


0.39 


-0.32 


0.37 


26.16 


3.37 


GSC 00154-00785 


06 33 36 


+05 16 41 


+ 1.87 


1.90 


-4.26 


2.60 


28.2 


1 


GSC 00154-01871 


06 33 43 


+04 55 54 


-1.16 


1.78 


-3.21 


1.70 







4 FIELD STARS OR CLUSTER [MEMBERS? 

4.1 Proper motion study and radial velocity 

To establish membersliip tlirough a proper motion study, NGC 2244 
is not an easy cluster to investigate, due to both its rather large 
distance and the similarity between the proper motions of stars 
in the cluster and those in the field, as is illustrated in the atlas 
of NGC 2244 published by iKharchenko et alj ( |2005| ). According 
to this atlas, the cluster member proper motions are concentrated 
around a mean of (- 1 .42, -0.32) mas/yr, compared to the mean for 
field stars in the same region of the sky of (+0.92, -6.55) mas/yr 
with a dispersion for the cluster stars of cTu ^^ = 4.45 ± 0. 15 mas/yr 
jKharchenko et al.l2004l : IChen et al.ll2007l) . 

Table [3] shows a comparison between the cluster mean values 
of proper motion and radial velocity for NGC 2244 and those of the 
two target stars: GSC 00154-00785 an d GSC 00154-018 71. From 
their independent proper motion study. Ichen et al. U2007h estimate 
membership probability for GSC 00154-00785 to be 50%, and for 
GSC 00154-01871, 84%. 

For GSC 00154-00785, the measurements of both fi„cos(d) 
and fis agree with those of the cluster within 2cr. The agreement 
of fiaCos(S) values improves if one consid ers only the 30 cluster 
members listed by Kharchenko et alj | |2004|) with measurement un- 
certainties smaller or comparable to those of GSC 00154-00785. 
In that sub-sample, the dispersion of /i„cos((5) is approximately 2.5 
mas/yr. The agreement between the cluster mean values and those 
of the two stars is as satisfactory for GSC 00154-00785 as it is for 
the most precisely measured stars in the cluster. 

The agreeme nt between the mean cluster radial velocity 
dKharcheiiko et alj [20 05) and the value given for GSC 00154- 
01871 bv lNiemczur^ foriv. commj) strengthens the case for mem- 



bership of this star. However, we do not include this star's ra- 
dial velocity in Table [3] because the star is a spectroscopic binary 
dHuang & Giesi2006l) and we do not have enough spectra over time 
to be able to provide a very accurate mean radial velocity of the 
binary system. For GSC 00154-01871, we arrive at the same con- 
clusion as for GSC 00154-00785, but even more strongly, since the 
star's value of fJaCos(6) agrees with the cluster mean to within only 
lo-. 



4.2 Reddening and CMD analysis 

In Fig.|6]we show the positions of the two stars relative to the cen- 
tral region of NGC 2244. 

Applying the calibration bv lCrawford jl978l) to the Stromgren 
colour indices gives E(b - >■) = 0.245 and M,, = -0.09 ± 0.4 
for GSC 00154-00785 and E{b - y) = 0.379, M,. = 0.15 ± 
0.4 for GSC 00154-01871. The considerably lower reddening of 
GSC 00154-00785 is understandable given its position in the outer 
parts of the region of NGC 2244 With V magnitudes of 10.93 and 
1 1 .30 respectively ( lHandlerl20"Tll) . we obtain distances of 980±400 
and 800 + 350 pc for GSC 00154-00785 and GSC 00154-01871, re- 
spectively. 

The distance o f NGC 2244 is estim ated to fall in the range 
of 1.4 - 1.7 kpc (see lBonatto & Bicall2009l and references therein). 
Our Stromgren results put the two stars at roughly half the cluster 
distance and thus argue against cluster membership. However, we 
point out that the distances we derive for the stars are within lo" of 
the lowest estimate of the cluster distance. 

Words of warning: If either star is surrounded b y circum- 
stellar matter, this would redden the stars, making the ICrawfor j 
( Il978l) calibration routine invalid. However, neither of the two stars 
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Figure 4. Upper panel: Region of interest of tlie amplitude spectram of 
GSC 00154-00785 (blue) with the identified frequencies (black). Lower 
panel: Amplitude spectrum in log - log space. Shown are the power-law 
fit to the data (green solid line) and the 3.8 o" limit (dashed line) 



shows emission lines in its spectrum /Niemczural Ipriv. commj) . 
GSC 00154-01871 is a spectroscopic binary and therefore its ab- 
solute magnitude derived from uvby photometry cannot be trusted. 
For B stars, the latter is calibrated via H/3 whose value is affected 
by the presence of a companion. 

The calibration routine implicitly assumes that the two stars 
are main-sequence (MS) stars. If these stars are in their pre-main 
sequence (PMS) stage, we argue that the Crawford (1978) routine 
can be safely applied because (a) there is no evidence for circum- 
stellar emission, and (b) the colour excesses from the uvby photom- 
etry are not anomalous. 

Figure |5] is a dereddened colour- magnitude dia g ram (CMP) 
contai ni ng stars from th e sampl es of Heisei jl977h ; Perez et al] 
(Il989h: iKaltcheva et al.l jl999l) ; iHensberge bOOd) and 

iHandleil ( 120111) . These have a membership probability > 85% 



(cros s -checked against Marshall et al. 1982; Kh archenko et al.l 



l2005l : iBaumgardt eTZI |2000| and IChen et al.. .20070 . The dered- 
dened co l our in dex {b - y)fi and Vo were obtained using the 
ICraw ford* ('1978") calibration as described above. The unsuspicious 
position of GSC 00154-00785 and GSC 00154-01871 in the CMD 
as well as their location in the field of NGC 2244 supports cluster 
membership of the two stars. 

4.3 Conclusion of this section 

Neither the cluster membership determinations nor the distance es- 
timates are iron-clad, and they suffer from qualifiers and caveats. 
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Figure 5. Zoom into the dereddened colour-magnitude diagram (CMD) 
for cluster members of NGC 2244. Positions of GSC 00154-00785 and 
GSC 00154-01871 are superimposed (red stars). 
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Figure 6. POSS-I image of NGC 2244 observed at a wavelength of 
~ 0.5/im. Green circles denote the positions of the two guide stars, 
GSC 00154-00785 and GSC 00154-01871. 



Although, membership of GSC 00154-00785 and GSC 00154- 
01871 is possible, we are unable to confirm it. GAIA parallaxes 
will settle this issue. 



5 MODELS 

5.1 Evolutionary models 

GSC 00154-00785 

The effective temperature and surface gravity derived in Sect. 13.21 
are plotted in the bottom panel of Fig.|7]with main-sequence (MS) 
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evolutionary track^ (solid lines) and pre-main-sequence (PMS) 
tracks (dashed lines). This diagram indicates that GSC 00154- 
00785 has a mass of 4.5 ± 0.5 no matter whether it is an MS 
or PMS star. 

We can estimate the age of the star by comparing its position 
on the log^ - log T^ff diagram with the evolutionary tracks. If the 
star is assumed to be a MS star, its age is 20 - 50 Myr. 

Assuming it to be a field star, thus using the absolute mag- 
nitude derived from Stromgren photomet ry (Sect. 14.2b . and ap- 
plying a bolometric correction of -1.14 jFloweilll996l) . we find 
M\ya\ = -1.23±0.4 which corresponds to a luminosity of log L/Lq = 
2.42 ±0.16. However, if GSC 00154-00785 is a cluster member, 
adoptin g a dist ance modulus of {m - M)o = 11.1 for NGC 2244 
(see e.g. lPark & Suna.200Z and references therein) we find Mboi = 
-2.27, corresponding to a luminosity of log L/Lq = 2.8 ± 0.25. 

GSC 00154-01871 

The position of GSC 00154-01871 in Fig.|7]is also based on the 
parameters derived in Sect. 13.21 The location of this star in Fig. |7] 
suggests that this star has a mass of 5.0 ± 0.5 M©, whether it is an 
MS or PMS star. If the star is assumed to be an MS star, its age is 
70 - 80 Myr. 

We estimate the luminosity as abo ve for GSC 0154-00785. 
With a bolometric correction of -1.02 CFlowej|l996) , we find for 
GSC 00I54-0187I the following values: Mb^, = -0.87 and a lumi- 
nosity of logL/Lo = 2.28 ± 0.16. On the other hand, GSC 00154- 
01871 has a stronger case for cluster membership. If we adopt the 
same distance modulus as for GSC 00154 -00785 above and after 
applying a bolometric correction of - 1 .02 l lFloweJl99d) we find a 
bolometric magnitude of Mboi = -2.45 and a corresponding lumi- 
nosity of log L/Lo = 2.91 ± 0.25. 

If we assume that GSC 00154-01871 is a field star, the evolu- 
tionary models that best fit this star in a log L - log Jeff plot con- 
flict with what is expected from the log g - log Teff diagram shown 
in Fig. [T] On the other hand, the stellar masses derived from the 
evolutionary models are consistent in both diagrams if the star is 
assumed to be a cluster member. Additionally, if the star is in the 
field, its luminosity would be only marginally consistent with mod- 
els of g-mode excitation. Both points are circumstantial evidence 
to support the membership of GSC 00154-01871 in NGC 2244. 

Implications: MS or PMS? 

If these stars are on the MS, the estimated ages are much larger 
than the age range fo r the open cluster NGC 2244 (0.2 - 6 Myr 
iBonatto & Bicall2009l) . This argues against a cluster membership 
if they are MS stars. However, if the stars are still in their PMS 
stage, the situation is different. The dashed lines in Fig. [7] show 
PMS evolutionary tracks. The initial model for each PMS track has 
been obtained by accreting mass starting from a I Mq model on 
the Hayashi track. As seen in Fig. [T] both stars can be identified 
as PMS stars with masses near 4.5 Mq. The ages are estimated for 
GSC 00154-01871 and GSC 00154-00785 to be about 0.2 Myr and 
0.4 Myr, respectively, after termination of mass accretion. If the 
stars are indeed cluster members, then their location in the HRD 
implies that they must be PMS stars. In the following, pulsational 
analyses, we consider the g-mode pulsation spectra in both MS and 

^ Evolutionary models were computed by a standard Henyey type code. 
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Figure 7. Estimated parameters for GSC 00154-00785 (filled symbols) and 
GSC 00154-01871 (open symbols) are compared with evolutionary models 
in the HR diagram (top panel) and in the log g - log T^ff diagram. The two 
positions for each star in the top panel correspond to the two kinds of esti- 
mates for the luminosity discussed in the text: circles denote the luminosity 
assuming the stars to be cluster members, triangles denote the field star as- 
sumption. Dashed lines are PMS evolution tracks, while solid lines are MS 
and post-MS evolution tracks for each stellar mass indicated (in solar units) 
along each track in the top panel. Dotted lines indicate the approximate 
locus of the birth line. 



PMS models. (A detailed discussion about g-modes in PMS stars 
can be found in Appendix lAl) 



5.2 Pulsation models 

Fig. [8] shows oscillation frequencies of excited low-degree I = \ 
and 1 = 2 modes along the evolutionary tracks of non-rotating 
4.5 Mq models in MS (left panel) and in PMS (right panel) stages. 
We have adopted a standard che mical composition {X = 0.7, 
Z = 0.02) with OP opacity tables jBadnell et aT]|2005l) . In a ca- 
sual inspection of Fig. [8] no significant differences in the range 
of excited g-mode frequencies are evident between MS and PMS 
models. Both types of models with similar effective temperatures 
excite g-modes with similar periods, because the excitation occurs 
in outer layers of the star. Before a convective core appears in the 
PMS stage (around the peak in luminosity), the amplitudes of g- 
modes in the deep interior are much larger than in MS models, 
although the number of nodes (a few tens) are not large enough 
to significantly dissipate g-mode oscillations as in the post-main- 
sequence models. The g-mode spectra for PMS models are simpler 
than those for MS models because of the absence of mode cross- 
ings associated with the gradient in mean molecular weight found 
after the ZAMS. 

The range of excited frequencies shifts downward as the ef- 
fective temperature decreases, and for a given effective tempera- 
ture the frequency range is slightly smaller for PMS models. Those 
trends are due mainly to the difference in stellar radii of the mod- 
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Figure 8. Frequencies of excited liigli-order, low-degree t = \ (blue triangles) and t = 2 (red squai'es) g-modes versus eff ective temperature for the 4.5 M© 
evolution models in the MS (left panel) and in the PMS (right panel) phases. Horizontal dotted and dashed lines indicate frequencies detected in GSC 00154- 
00785 and GSC 00154-01871, respectively. The cross and the horizontal Hne in the lower part of each diagram indicates the estimated effective temperatures 
with error bars for GSC 00154-00785 and GSC 00154-01871. 



els; the frequency of a g-mode is lower for a star of larger radius. 
Fig. [8] shows that, for both MS and PMS models, the frequency 
range of excited £ = \ and 2 modes covers well the ranges of os- 
cillation frequencies detected in GSC 00154-01871 (dashed lines) 
and GSC 00154-00785 (dotted lines). 

If GSC 00154-00785 and/or GSC 00154-01871 are PMS stars, 
evolutionary period (frequency) change is expected to be more 
rapid than in the MS stage. Fig.|9]shows the rates of period change 
dPjdt for selected g-modes in the 4.5 Mq model during PMS evo- 
lution. 

At the effective temperature of GSC 00154-00785 (log T^g = 
4.164), the periods are expected to increase with time, while at the 
effective temperature of GSC 00154-01871 (log T^g = 4.143), they 
decrease. The period decrease is caused by the PMS contraction, 
while the increase is caused by the fact that Brunt-Vaisala fre- 
quency decreases deep in the star with the development of a con- 
vective core which first appears at log Teff - 4.15. 

If GSC 00154-00785 and GSC 00154-01871 are PMS stars, 

Fig.|9]indicates that the period change rates is ~ 1 sec/yr and 0.2 

sec/yr, respectively. Such period changes might be detectable by an 
(O-C) diagram based on long baseline observations. 




4.16 4.14 
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6 CONCLUSIONS 

MOST has revealed the SPB nature of stars (GSC 00154-00785 
and GSC 00154-01871) in the same part of the sky as the very 
young open cluster NGC 2244. The main pulsation periods of these 
multi-periodic SPB stars are near 2 days. The ranges of pulsation 
frequencies of the stars is consistent with nonradial g-modes low 
degree (£ = i and C = 2) excited by the kappa mechanism at the Fe 
bump. 

The effective temperatures and surface gravities we derive 



Figure 9. The rates of period change for some g-modes in the 4.5 Mq model 
during the PMS evolution. Plotted modes are gi^, g2o, gis ™d g30 for t = I 
and 2. (Longer period modes tend to have larger \dP/dt\.) Filled triangles 
and squares indicate modes for C = I and 2, respectively to be excited. 

from Stromgren indices and spectroscopy indicate the two stars 
have masses of about 4.5 - 5.0 Mq. If GSC 00154-00785 and 
GSC 00154-01871 are MS stars, our estimates for their ages are 
20 - 50 Myr and 70 - 80 Myr, respectively, much older than the 
cluster NGC 2244 (0.2 - 6 Myr). If the two stars are still in the PMS 
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stage, they would be a few tenths of Myr after the end of the pro- 
tostellar accretion phase. Therefore, if these stars are proven mem- 
bers of NGC 2244, then almost certainly they are PMS stars. Ex- 
isting proper motion and radial velocity properties of the two pul- 
sators cannot establish reliably their membership in the NGC 2244 
cluster. GAIA parallaxes and proper motions should eventually re- 
move this uncertainty. 

If the membership to NGC 2244 can be indeed conflrmed, 
GSC 00154-00785 and GSC 00154-01871 would be the first known 
PMS SPB pulsators. This possibility has alerted us to begin a more 
concerted search with MOST and other space photometry missions 
for other SPB pulsators among PMS cluster and field stars. 
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APPENDIX A: INSTABILITY RANGE FOR G-MODES 
IN PMS STARS 

Since PMS evolutionary tracks of intermediate masses (3 <, 
M/Mq < 6) enter into a g-mode instability region (i.e., the 
SPB instability region) caused by the Fe opacity bump al T x: 
1.5 X 10^ K, certain g-modes can theoretically be excited in some 
intermediate PMS stars. We have examined the instabilities of 
nonradial g-mode oscillations of ^ = 1 and 2 for PMS models 
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Figure Al. Models in which g-modes are excited are indicated by sym- 
bols (A and □ for C = \ and 2 modes, respectively for PMS models, and 
X for both Cs for main-sequence models ) along evolutionary tracks with 
various masses. Red solid lines indicate PMS models, while blue dashed 
lines main-sequence and post-main sequence models. Dotted line indicates 
evolutionary tracks during mass-accreting proto-stellar evolution. 

with masses ranging across 2.7 M© to 6 Mq. The initial models 
for the evolution of each mass were obtained by accreting mass 
at a rate of 10~^ Mgyr"' starting from a fully convective 1 Mq 
model. Those initial mo dels are roughly located on the birth-line 
of PMS stars l lPalla & St abler 1993). 

The stability results are shown in Fig. lAll where red trian- 
gles and squares indicate the positions of models in which g- 
modes are excited. Red lines indicate the evolutionary tracks of 
PMS stars, while blue lines are main-sequence and post-main- 
sequence tracks. Main-sequence models unstable to g-modes 
(corresponding to what would be currently recognised SPB stars) 
are highlighted by blue crosses. 

Because the Brunt-Vaisala frequency becomes very high in the 
core of a post-main-sequence star, g-modes are damped due 
to dissipation unless a shell convection zone appears, as in 
massive stars iSaio et aL I l200d lGautschvl[2009h . This is why 
the instability range for prototype SPB stars is confined to the 
MS band (Gautschv & Saio 1993; Dziembowski et al. 1993). In 
PMS stars, however, g-modes survive even in the phase with- 
out core convection, because the Brunt-Vailala frequency in the 
interior is not as high as in a post-main-sequence star. The spa- 
tial wavelength of the eigenfunction is therefore too long for the 
dissipation to be important. (The number of nodes of excited g- 
modes is well below 100.) As a result, the g-mode instability 
range for PMS stars (red symbols in Fig. lAlt is broader than that 
for MS SPB stars (blue symbols). The models predict PMS SPB 
stars which overlap with and extend beyond the main sequence 
band in the HR Diagram. 



